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1085,  United  States 

Abstract 

The  present  paper  addresses  estimation  of  global  deformation  of  a  continuous-fiber  brittle 
matrix  composite  over  a  length  much  larger  than  the  largest  microstructural  dimension  of 
the  composite,  such  as  the  fiber  diameter  and  the  fiber  spacing.  Initial  microcracks  and 
progressive  microcracking  of  the  matrix  and  interfaces  under  time-dependent  load  and 
temperature  are  considered.  The  model  also  includes  consideration  of  inelastic 
deformation  of  one  of  the  constituents.  In  this  paper  the  equations  of  the  model  are 
developed  and  applied  to  model  the  stress-strain  behavior  of  several  SiC/SiC  CMCs  and 
the  creep  behavior  of  an  oxide/oxide  CMC.  Comparison  with  test  data  shows  that  the 
model  is  able  to  capture  the  wide  range  of  deformation  behavior  seen  in  these  CMCs. 

Introduction 

A  distinguishing  characteristic  of  the  class  of  ceramic-matrix  composites  is  that  the  stress 
fa)  and  time  ft)  dependent  strain  (s/a,  t}),  is  partly  due  to  micromechanical  damage  in 
the  form  of  matrix  cracks  and/or  interfacial  separations  and  partly  due  to  inelastic 
deformation  of  a  constituent.  Detailed  observations  on  such  composites  have  been 
presented  by  Karandikar  and  Chou  [1],  Zhu,  et.  al.  [2]  and  others.  Time  dependence  of 
damage  and  deformation  mechanisms  due  to  creep  and  due  to  chemical  reaction  of 
constituents  with  surrounding  environment  have  been  discussed  by  DiCarlo  [3], 
Morscher  [4],  Eckel  [5]  and  others.  Based  on  observations  reported  in  the  literature, 
Figure  1  is  a  schematic  representation  of  microscopic  damage  and  deformation 
mechanisms  and  their  manifestation  on  macroscopic  strain,  s/a,  t},  of  the  composite. 

Under  monotonically  increasing  load  over  a  short  time,  general  tensile  stress-strain 
behavior  of  the  composite  is  shown  in  Figure  1(a).  It  shows  an  approximately  linear 
portion  (regions  1  and  2),  followed  by  a  gradually  decreasing  slope  ( da/ds )  in  region  3, 
and  subsequent  stiffening  in  region  4.  Relative  size  of  each  region  varies  with  material 
system.  For  example,  in  some  materials,  composite  rupture  occurs  without  any 
perceptible  stiffening.  In  others,  presence  of  microcracks  in  as-fabricated  condition  can 
result  in  virtual  absence  of  regions  1  and  2,  resulting  in  a  stress-strain  curve  that  is 
nonlinear  throughout. 

The  stress  at  the  end  of  region  1  is  the  ‘proportional  limit’,  where  a  discemable  deviation 
from  the  initial  linear  behavior  occurs.  Up  to  a  certain  stress,  usually  less  than  the 
proportional  limit,  the  composite  suffers  no  additional  damage.  However,  pre-existing 
defects  and  damage  (such  as  matrix  porosity,  interfacial  disbonds  and  microcracks)  may 
be  present,  affecting  the  initial  slope  as  well  as  the  proportional  limit. 
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Transverse  microcracks  first  appear,  usually  within  the  transverse  tows  in  woven 
composites,  at  a  stress  less  than  the  proportional  limit  [1],  shown  as  start  of  region  2  in 
Figure  1(a).  Some  microcracks  may  be  “gaps”  among  fibers  within  transverse  tows  due 
to  insufficient  matrix  material  penetration  into  the  tows.  Such  gaps  can  open  up  under 
very  low  load.  The  onset  of  microcracking  may  also  be  influenced  by  composite 
processing  induced  residual  stresses  due  to  mismatch  of  coefficient  of  thermal  expansion 
(CTE)  among  the  constituents.  Under  increasing  load  microcracks  enlarge,  additional 
microcracks  form,  and  splits  (longitudinal  cracks  or  interfacial  disbonds)  appear  in  region 
2  as  the  applied  stress  approaches  often  a  distinct  knee  in  the  stress-strain  curve.  Higher 
stress  (region  3)  causes  accelerated  microcracking  accompanied  by  enlargement  of  splits. 
Some  longitudinal  fibers  may  fracture  toward  the  end  of  region  3  as  increasing  matrix 
damage  causes  accelerated  load  transfer  to  the  fibers.  In  region  4,  applied  load  is  carried 
largely  by  the  longitudinal  fibers,  which  eventually  fracture,  usually  in  rapid  succession, 
resulting  in  composite  rupture. 

As  shown  in  Figure  1(a),  with  insufficient  time  for  significant  time  dependant 
deformations  to  have  occurred,  unloading  from  various  stress  levels  and  reloading  is 
approximately  linear,  with  little  or  no  pennanent  deformation.  However,  hysteresis  loops 
occur  in  some  composites,  especially  at  elevated  temperatures. 

Figure  1(b)  depicts  strain-time  response  of  the  composite  under  a  constant  applied  load 
and  at  a  sufficiently  high  temperature  for  constituent  creep  and/or  chemical  reaction  of 
constituents  with  surrounding  environment  to  be  significant.  The  sequence  of 
micromechanical  damage  progression  with  time  is  conjectured  to  be  similar  to  that  under 
increasing  tensile  load,  but  observations  reported  in  the  literature  are  less  definitive 
because  of  difficulties  associated  with  experiments  at  high  temperatures.  At  constant 
applied  load,  strain  increase  can  be  due  to  creep  of  one  or  more  constituents,  interfacial 
sliding,  additional  microcracks  and  crack  growth  in  the  matrix,  and  discrete  fiber 
fractures  induced  by  reaction  with  environment.  For  example,  in  some  composites, 
oxidation  of  the  matrix  can  lower  its  strength,  resulting  in  gradually  increasing  number  of 
microcracks  under  constant  applied  load. 

Depending  on  the  composite,  a  portion  of  the  time  dependent  strain  may  be  permanent. 
Creep  of  one  or  more  constituents,  irreversible  sliding  at  interfaces,  and  “wedging”  of 
microcracks  by  reaction  products  are  among  plausible  causes  for  permanent  defonnation. 

There  are  several  reported  approaches  in  the  literature  for  modeling  micromechanical 
damage  and  defonnation  mechanisms  in  the  types  of  composites  presently  considered. 
The  model  of  Aveston,  Cooper  and  Kelly  [6]  has  served  as  the  basis  of  a  number  of 
approaches,  with  significant  improvements  and  generalizations  by  Marshall,  Cox  and 
Evans  [7],  Curtin  [8]  and  others.  A  detailed  mechanics  analysis  of  matrix  cracks  by 
Budiansky,  Hutchinson  and  Evans  [9]  and  of  fiber  creep  by  Begley,  Evans  and 
McMeeking  [10]  have  further  enhanced  the  understanding  of  damage  and  deformation  of 
brittle  matrix  composites.  More  recently,  the  stress-strain  behavior  of  CMCs  has  been 
modeled  using  a  phenomenological  approach  [11]  and  using  continuum  damage-based 
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approach  [12-15].  None  of  these  approaches  however  considers  inelastic  deformation  of 
the  constituents,  and  in  most  instances,  the  models  require  the  testing  of  a  range  of 
different  laminate  configurations  in  order  to  obtain  the  inputs  required  for  the  application 
of  the  model. 

The  approach  in  the  present  work  builds  on  several  concepts  and  ideas  found  in 
references  [6-10],  and  introduces  extensions  and  enhancements  aimed  at  improved 
modeling  of  deformation  and  damage  mechanisms  in  CMCs  relevant  to  current 
engineering  applications.  Special  attention  is  given  to  developing  models  that  include 
consideration  of  at  least  all  the  dominant  damage  and  deformation  mechanisms,  while 
keeping  the  model  reasonably  simple,  computationally  efficient,  and  minimizing  the 
number  of  laboratory  tests  needed  to  generate  input  for  the  model. 

Model  Derivation 

To  simplify  presentation  of  the  model,  let  us  consider  a  two-phase  composite,  with  fiber 
and  matrix  as  the  only  constituent  materials.  In  such  a  material  the  longitudinal  and 
transverse  stiffness  of  the  composite  can  be  approximated  by  the  rule-of-mixtures  (ROM) 
and  the  inverse  ROM  respectively  [16].  Furthermore,  let  us  prescribe  the  matrix  material 
to  be  linear-elastic.  Let  us  express  the  inelastic  strain  rate  of  the  fiber  material  by  the 
following  three-dimensional  (k  and  /  =  1,3)  rate  form  of  the  Prandtl-Reuss  equations: 

& kl(  f  )  =  ^kl(  f  )■>  (1) 

where  X  is  a  scalar  function  of  state  variables.  If  the  fiber  is  in  pure  axial  tension,  with  all 
other  stress  components  being  zero, 

4=|  A(Tf  (2) 

Both  the  matrix  and  the  fiber  materials  are  assumed  homogeneous.  It  is  implied  that  any 
porosity  in  the  matrix  is  uniformly  distributed  and  the  largest  pore  size  is  small  compared 
to  the  woven-architecture  dimensions  such  as  tow  diameters  and  tow  spacing.  Therefore, 
porous  matrix  can  be  approximated  as  homogeneous,  albeit  with  mechanical  and  thermal 
properties  different  than  a  fully  dense  matrix  material.  For  example,  the  Young’s 
Modulus  of  the  matrix  can  be  estimated  as:  EM  =  (1  -  p )EM(dense) ,  where  p  is  the  total 

volume  of  pores  divided  by  the  volume  of  the  matrix  in  the  composite.  Similarly,  other 
mechanical  and  thermal  properties  of  a  porous  matrix  can  be  estimated.  Of  course,  more 
accurate  models  to  account  for  porosity  are  possible.  But,  for  the  present,  these  simple 
estimates  suffice. 

Composite  Without  Damage 

Let  us  first  consider  a  strip  of  length  2L,  width  IV u,  and  uniform  thickness  d,  of  a 
composite  without  damage,  and  having  v/  volume  of  fibers  per  unit  volume  of  the 
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composite  (i.e.,  the  average  fiber  volume  fraction  is  vj).  The  field  quantities  specific  to 
the  undamaged  composite  are  identified  with  the  subscript  “U”.  The  length  is  along  one 
of  the  fiber  directions.  Each  edge  of  the  strip  is  displaced  an  amount  (8p)  by  a  force  P\j 
acting  normal  to  the  cross-section  (Wu'd).  Thus,  the  average  stress  on  the  composite  is 
( 7u=Pv/(W-d ')  and  the  total  strain  is  Su=8u/L.  Each  dimension  of  the  composite  strip  is 
large  compared  to  microstructural  dimensions  of  the  composite.  In  the  case  of  a  woven 
composite,  microstructural  dimensions  include  dimension  of  fiber-tow  and  the  amplitude 
and  wavelength  characteristics  of  the  weave. 

If  the  composite  is  unidirectional,  the  applied  stress  <Ju  is  supported  by  the  matrix  and  the 
fibers.  In  a  woven  composite  with  orthogonal  tows,  <Ju  is  shared  by  the  longitudinal 
fibers,  the  matrix,  and  the  transverse  fibers.  If  the  wave-length  to  amplitude  ratio  of  the 
weave  is  large,  such  a  woven  composite  can  be  viewed  as  a  balanced  cross-ply  laminate, 
i.e.  a  (0/90)s  laminate,  with  equal  volume  of  fibers  (v//2)  in  the  load  (0°)  direction  and 
the  transverse  (90°)  direction.  In  that  case,  the  woven  composite  can  be  approximated  as 
a  unidirectional  composite  having  Vf/2  fiber  volume  fraction,  and  a  matrix  composed  of 
the  actual  matrix  material  (M)  and  the  transverse  fibers.  Let  us  denote  this  ‘effective’ 
matrix  by  m,  and  half  the  fiber  volume  fraction  by  Vf,  such  that,  for  a  unidirectional 
composite,  m  is  the  same  as  M  and  Vf=Vf. 

In  woven  composites,  fiber-tows  may  contain  matrix  (M),  voids,  and  individual  fibers, 
with  a  different  fiber  volume  fraction  (vfw)  within  the  tows  than  the  average  volume 
fraction  (vp)  in  the  composite.  In  that  case,  v/=vw.y/ -w,  where  vw  is  the  volume  of  tows  per 
unit  volume  of  the  composite. 

Assume  that  under  time-dependent  applied  stress,  <jjj  (t),  the  fibers  (f)  and  the  effective 
matrix  (m)  undergo  the  same  total  average  strain  rate  ( sv  ),  i.e., 


£U  ~  £Uf  ~  £Um 


(3) 


At  strain  (£v)  equal  to  zero,  the  fiber  and  the  matrix  may  have  initial  (residual)  stresses 

Rf  and  Rm,  respectively,  induced  by  the  composite  fabrication  process.  For  example,  the 
residual  stresses  can  be  due  to  unequal  thennal  contraction  characteristic  of  the  two 
materials  during  cool-down  from  some  stress-free  temperature  ( Tp)  to  a  reference 
temperature  (Tr).  Simple  estimates  of  such  thennal  residual  stresses  can  be  obtained 
using,  for  example,  thermoelastic  analysis  results  of  Gatewood  [17]  for  concentric 
cylinders  of  dissimilar  materials. 

The  strain  rate  in  the  fiber  consists  of  elastic,  thennal  and  inelastic  strain  rates,  i.e., 


£u  ~  D, 


auf  Rf 
V  Ef 


+  af(T-Tr) 


2 

+  j  Ajf&uf  ’ 


(4) 
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where,  Dt  denotes  derivative  of  the  quantity  within  parenthesis  with  respect  to  time,  and 
T=T{t},  Ef=Ef[T}  and  af^atfT}  are  the  temperature,  longitudinal  fiber  modulus  and  fiber 
CTE  respectively. 


Eq.  (4)  can  be  rearrange  to  give: 


^uf  ~  Ef 


£v  -  cTufE>,(l / E /)- E>t  af(T-Tr)-^f 

{  EfJ 


(5) 


The  matrix  being  linear  elastic, 


&U  —  Dt 


T  _  _  n  \ 

- -  +  a(T-Tr ) 


v 


j 


(6) 


Here  Em  and  am  are  the  Young’s  modulus  and  CTE  of  the  effective  matrix  respectively. 
As  mentioned  earlier  in  a  unidirectional  composite  system  the  effective  matrix  is  the 
same  as  the  actual  matrix,  but  in  a  balanced  cross-ply  system  the  Young’s  modulus  and 
CTE  of  the  effective  matrix  have  to  be  determined.  Expressions  for  these  quantities  in 
terms  of  the  longitudinal  (EJ)  and  transverse  (Ejf)  modulus  of  the  fiber,  Poisson’s  ratio 
(d/)  and  CTE  (aj)  of  the  fiber,  the  modulus  (EM),  Poisson’s  ratio  (%)  and  CTE  (aM)  of  the 
actual  matrix  and  the  volume  fraction  of  the  fiber  (v/)  are  given  in  the  appendix. 

Note  that  with  the  assumptions  already  cited, 

O'u  =  VfaUf  +  V,naUm  »  (7) 


Eu=VfEf  +  VmEm 


(8) 


and,  the  matrix  being  linear  elastic, 

°Um  =  E,n (£U  -  am (T~Tr))+  Rm  •  (9) 


Note  that,  depending  on  how  the  volume  of  voids  in  the  matrix  and  the  tows  is  accounted 
for,  the  volume  fraction  of  the  matrix  (Vm)  may  not  equal  (1-Vf),  and  that  f/  and  Vm  can  be 
functions  of  time  and  temperature. 


Therefore, 


°Uf  ~ 


_  au 


1-jEjt:,  ajT  /  ))•  AfJ 
Vf 


(10) 


Substituting  Eqs.  (5)  and  (1 1)  in  Eq.  (6),  we  get, 
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%  = 


CTv-[E,(su-aJT-Tr))+R,]v„ 

V, 


^  C ' f\if  +  CT 


+  c,  +C  D, 

fin  m  t 


o, 


\V  E  , 

\  m  my 


(11) 


In  the  above  equation. 


Cf  =  VfEf/Ea ’ 

C,n  =  V,nE,nlEU  ’ 

r 


=  CfD, 

\af(T-Tr)- 

fM 

1  T7 

\ 

+  CmDt 

iam(T-Tr)- 

y)) 

\  F 

V 

fEf  ), 

/ 

V 

\^m  J  J 

and 

CT  =  CfDt 


\EfJ 


CD, 


K) 

V  E  , 

\  m  my 


Eq.  (11),  together  with  an  appropriate  initial  condition,  can  be  solved  for  sv{t, 

For  example,  one  can  choose  the  initial  time  t=to,  temperature  T=T0,  and  applied 
stress  cXy  =  cxu0  such  that  the  corresponding  strain  is  purely  elastic  and  is  given  by: 


on 


'U  I T=Tn 


(12) 


Except  for  some  special  cases,  the  above  initial  value  problem  requires  a  numerical 
solution  method.  One  case  amenable  to  exact  solution  is  when  the  composite  is  subjected 
to  a  constant  applied  load  at  constant  temperature  T=Tr ,  and  the  fiber  material  undergoes 
inelastic  (creep)  in  uniaxial  tension  according  to  the  following  relation: 


s'f  =  Aonf 


(13) 
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Then,  Ayf  =  —Ao’^1  and,  for  the  initial  condition  given  by  Eq.  (12),  the  solution  is: 


%  = 


Vf 


cymEm 


f  ^ 

( 

C i 

^E  +  Rf 

_ 

vf  f  J 

\  f  J 
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where,  C,  =  (a  Cf )/  "  and  C2  = 


C, 


Aji-cj  ' 


-\i- 


v  F/ 


-  +  R, 


Composite  With  Fully  Cracked  Matrix 

Let  us  now  consider  a  unit  cell  of  length  2L,  width  Wb,  and  uniform  thickness  d,  of  a 
composite  in  which  all  matrix  has  cracked.  The  field  quantities  specific  to  this  damaged 
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composite  are  identified  with  the  subscript  “B”.  The  composite  is  subjected  to  loading 
similar  to  that  considered  previously  for  undamaged  composite. 

It  is  assumed  that  once  a  crack  in  the  matrix  nucleates,  it  rapidly  propagates,  in  directions 
nonnal  to  the  load  direction,  to  cover  the  entire  width  and  depth  of  the  matrix,  arresting  at 
or  near  the  fiber-matrix  interface,  and  causing  a  split  of  length  a  along  the  load  direction 
(Figure  2).  In  line  with  the  split,  there  is  a  distance  s  over  which  shear  stress  in  non-zero. 
The  effective  Young’s  Modulus  (E/12)  of  region  B2  of  length  (a+s)  is  VjEf. 

Strain-rate  of  Type  B  cell  can  be  expressed  as: 

£B  =  P£B2  +  (!  ~  P)£  B1+P{£B2  -  £B1 )  ( 1 5) 

in  which,  p  =  (a  +  s)/L.  If  the  matrix  is  brittle  or  if  there  fibers  within  a  tow  are  weakly 
bonded,  s  is  expected  to  be  much  smaller  than  a.  For  the  present,  let  us  assume  that  this 
is  the  case  and,  therefore,  p~  a/ L  . 


Note  that  region  B1  is  structurally  the  same  as  the  undamaged  composite  (U)  considered 
previously.  Therefore,  strain  rate  ( sB, )  of  this  region  is  given  by  Eq.(ll)  with  the 
subscript  UU”  is  replaced  by  UB1  ”,  and  assigning  crBI  =  <jb  . 


£B1 


\Em(£B,-<Xn,(T  ~Tr))+  KYn, 


2 

—Cf?iBif  +  CT 


- 

(  _  h 

+  C,  +C  D, 

fm  m  t 

®  B 

VE 

\  m  my 

(16) 


Strain  rate  in  region  B2  is  simply  the  strain  rate  in  fiber  under  stress  equal  to  <jB/Vf ,  i.e., 


'B2 


=  D,( 


^  +  af(T-Tr))  +  EB!ff- 


(17) 


Thus,  an  expression  forfy  is  found  by  substituting  Eq.  (16)  and  Eq.  (17)  in  Eq.  (15).  An 
expression  for  growth  rate  of  non-dimensionalized  split  length,  p,  can  be  assigned  based 
on  a  suitable  criterion.  Eq.  (15),  together  with  appropriate  initial  conditions,  can  be 
solved  for  sH  \t,  aB ,  T)  . 

Composite  With  Progressive  Matrix  Cracks 

We  now  consider  a  composite  strip  of  length  2L,  width  W,  and  uniform  thickness  d,  of 
which  any  cross-section  normal  to  the  load  contains  cells  of  equal  width,  some 
undamaged  (Type  U)  and  some  damaged  (Type  B)  as  shown  in  Figure  3.  Let  the  ratio  of 
the  number  of  Type  B  cells  to  the  total  number  of  cells  be  77 .  Let  the  strip  be  subjected  to 
applied  stress  a  =  P/(Wd )  such  that  su  =sB.  Then,  stress  acting  on  region  B  is  given 
by: 
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&b  =  \?-(l-V)<ru]lv 
and 


i)&u -fj 


f  a  -a,  ^ 


(18) 


(19) 


The  total  strain  and  rates  of  strain  for  regions  U,  B  and  the  composite  are  the  same,  i.e., 

£  =  £u=£b  (20) 

e  =  eu=eB  (21) 

Rearranging  Eq.  (11),  and  using  Eqs.  (20)  and  (21),  we  have, 


a,  = 


V  E 

m  m 


C_ 


£  ~ 


^u-[Em{s-aJT-Tr))  +  Rmym 

V, 


j  f  f^Vf  CT 


C fin  ' 


v  K„Em  j 


(22) 


Substituting  Eq.  (22)  in  Eq.  (19),  we  get, 
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(23) 


Substituting  Eqs.  (19)  and  (23)  in  Eqs.  (16)  and  (17),  we  get  expressions  for  sBI  and  sB2  . 
These  expressions,  when  substituted  in  Eq.  (15),  together  with  Eq.  (20),  provide  the 
following  expression  for  strain  rate  in  composite  with  progressive  matrix  cracks: 

£  =P^b2+(1-PMbi+P^b2-^bi)  (24) 


The  expression  for  e  includes  integrated  values  cr ,  T,  gv  ,  £\, ,  elB1  and  slB2 .  Therefore, 
the  solution  of  Eq.  (24)  requires  prescription  of  appropriate  initial  conditions  for  these 
variables.  Additionally,  expressions  for  p  and  // ,  and  corresponding  initial  conditions 
for  p  and  //  are  needed. 

Choice  of  p  and  77  Functions 

At  a  given  load,  the  damage  parameter  77  represents  the  fraction  of  the  total  number  of 
matrix  regions  between  fibers,  normal  to  the  load  direction,  that  have  cracks.  Thus,  using 
Weibull  statistics,  we  choose  the  following  function  for  77: 
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rj 


=  7o  +(l-*7<J  1-e 


-0.693 


=  0 


for  cr,„  >  0 


for  cr,„  <  0 


(25) 


In  the  above,  <Jm  is  the  maximum  of  the  stress  in  the  matrix  in  regions  XJ  and  Bl,  r/N  is  the 
Weibull  modulus,  r/0  is  initial  value  of  r/  and  cr0  is  median  strength  of  the  matrix.  At 
a m  =  a 0 ,  rj=0.5.  Both  qN  and  <j0  may  be  temperature  and  time  dependent.  For  example, 
<j0  may  decrease  with  time  due  to  the  environment-induced  degradation  of  the  matrix,  or 
with  number  of  load  cycles  due  to  fatigue. 

An  expression  for  the  second  damage  parameter  p  should  be  based  on  fracture  mechanics 
concepts.  In  so  far  as  the  split  is  induced  by  the  arrest  of  a  transverse  crack  in  the  matrix, 
it  can  be  argued  that  its  initial  length  is  proportional  to  the  stress  intensity  factors 
associated  with  the  matrix  crack.  Assuming  that,  on  the  average,  each  matrix  crack  is  of 
equal  length  -  say,  the  average  fiber  spacing,  then  the  associated  stress  intensity  factors 
are  linear  functions  of  only  um .  Thus,  it  can  be  argued  that  the  initial  split  length  is  a 

function  of  only  <jm .  However,  subsequent  growth  of  the  split  would  be  dictated  by  the 
stress  intensity  factors  associated  with  the  tip  of  the  split  itself. 

Even  if  it  is  possible  to  estimate  stress  intensity  factors  associated  with  a  growing  split, 
eventually  an  empirical  approach  would  be  needed  to  develop  growth  rate  criteria  for  the 
splits,  requiring  measurements  of  split  length  as  it  would  grow  under  load.  At  present, 
such  measurements  are  infeasible  even  under  ambient  conditions.  What  may  be  feasible 
is  to  infer  split  growth  rate  criteria  from  global  defonnation  characteristics  of  the 
composite. 

Considering  the  above  difficulties,  let  us  consider  an  opportunistic  approach  and  assume 
that  the  functional  form  of p  is  the  same  as  of  rj .  Specifically,  we  assume  that: 


P  =  Po  +(!-A)) 


1  -  e 


-Pc 


for  rj  >  0 


for  tj  <  0 


(26) 


In  the  above,  pc  and  pN  are  time  and  temperature  dependent  parameters  and  p0  is  initial 
value  of  p .  It  is  shown  later  that  the  above  choice  of p  allows  modeling  of  a  broad 
spectrum  of  global  defonnation  characteristics  of  composites.  By  interpreting  the  77  and 
p  parameters  as  measures  of  the  matrix  cracks  parallel  to  and  perpendicular  to  the  fibers 
respectively,  the  product  of  77  and  p  can  be  interpreted  as  “crack  density”  in  a  two- 
dimensional  sense. 
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In  summary,  the  determination  of  damage  state,  defined  by  p  and  77 ,  requires  inferring 
the  values  of  an,  r/N,  pc  and  pN  for  the  composite.  In  principle  it  should  be  possible  to 
find  (J0  and  >]N  using  statistical  distribution  of  strength  of  the  matrix;  at  least  in  cases, 
such  as  some  unidirectional  composites,  where  matrix  represents  a  single  monolithic 
material.  However,  in  general,  this  is  not  feasible.  Instead,  like  pc  and  pN  ,  <j„  and  ?;¥ 
must  be  inferred  using  a  composite’s  global  defonnation  characteristics. 

Progressive  Failure  of  the  Fibers 

As  will  be  shown  later  the  effect  of  matrix  cracking  as  modeled  by  the  77  and  p  functions 
of  equations  (25)  and  (26)  is  to  reduce  the  modulus  of  the  composite  and  thus  introduce 
nonlinearity  in  the  composite  stress-strain  behavior.  As  more  of  the  matrix  cracks  the 
composite  strain  increases  until  finally  the  fiber  strain  approaches  its  critical  value  at 
failure.  The  failure  of  the  reinforcements  in  turn  leads  to  the  ultimate  failure  of  the 
composite.  Thus  in  theory  it  should  be  possible  to  model  the  final  failure  of  the 
composite  using  the  failure  statistics  of  in-situ  fibers.  Using  a  Weibull  form  to  model  the 
failure  of  the  in-situ  fibers  we  can  write  the  following  expression  for  the  fraction  of 
unbroken  fibers: 


1//  =  exp 


Vc 


S~£n 


\£f  £o  J 


for  £  >  s0 


(27) 


where  £0  is  the  minimum  value  of  strain  where  the  fiber  starts  to  fail  and  y/c  and  y/n  are 
parameters  defining  the  Weibull  distribution.  As  before  by  setting  y/c  =-0.693,  £f 

becomes  the  median  failure  strain  of  the  in-situ  fibers.  In  the  absence  of  fiber  failure 
strain  data  we  can  estimate  the  remaining  Weibull  parameters  in  equation  27  from  the 
composite  failure  data.  The  ^parameter  is  then  used  to  progressively  reduce  the  volume 
fraction  of  the  fibers  in  the  composite. 

Solution  of  the  CMC  Model 

Exact  solution  of  the  CMC  model  described  by  equation  (24)  can  be  obtained  for  the 
special  case  of  a  two-phase  composite  in  which  the  fiber  and  the  matrix  materials  are 
both  linear  elastic.  For  this  case  non-linear  stress-strain  behavior  of  the  damaged 
composite  can  be  described  by  the  relation 


ct  =  E£  (28) 

where  E  is  the  effective  modulus  of  the  composite  containing  matrix  cracks  and  is  given 
by 
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E  =  (1  -  r/)Ec  +  t ! 


VfEfEc 


{l-p)VfEf+pE, 


(29) 


Here  Ec  is  the  modulus  of  the  undamaged  composite,  which  is  also  the  modulus  of  cell-U, 
given  by  equation  (8).  Equation  (28)  can  be  evaluated  in  closed-fonn  after  substituting 
the  expressions  (25)  and  (26)  for  the  rj  and  p  functions  respectively  and  by  noting  that  the 
matrix  stress  in  these  expressions  is  simply 

°m=EJs-aJT-Tr))  +  Rm  (30) 

where  Em,  am  and  Rm  are  respectively  the  Young’s  modulus,  CTE  and  the  residual  stress 
of  the  effective  matrix. 

A  more  general  application  of  the  foregoing  model  requires  numerical  solution  of  the 
equation  (24)  together  with  appropriate  initial  conditions.  In  the  present  work  the 
numerical  solution  method  used  is  the  fourth-order  Runge-Kutta  method  [18]. 

Figure  4(a)  shows  the  stress-strain  response  of  a  CMC  determined  using  the  finite- 
difference  numerical  solution  method.  The  solution,  obtained  using  two  different  values 
of  time-step,  At,  is  compared  with  the  exact  solution  detennined  using  equation  (3.34). 
The  values  of  the  material  constants  and  damage  parameters  used  to  determine  the 
solution  are  also  included  in  the  figure.  The  results  indicate  that  as  the  time  step  is 
decreased  the  numerical  solution  converges  to  the  closed-fonn  solution.  Similar 
comparison  of  the  damage  functions,  rj  and  p,  obtained  using  the  two  solution  techniques 
are  shown  in  Figure  4(b).  The  results  demonstrate  the  accuracy  of  the  finite-difference 
numerical  solution  method  used  to  solve  the  nonlinear  equations. 

Results  and  Discussion 

Figures  5,  6  and  7  show  the  room-temperature  (RT)  tensile  stress-strain  curves  of  three 
ceramic-matrix  composites  generated  using  the  present  model  and  the  corresponding 
laboratory  test  measurements.  Note  that  the  numerical  solutions  are  not  predictions. 
Rather,  with  the  exception  of  fiber  properties  and  architecture  specifics  of  the  composite, 
all  other  input  data  for  the  model  were  found  using  the  measurements  shown  in  the 
figures.  While  the  solutions  are  not  direct  curve-fits  of  the  measured  data,  they  are 
nevertheless  “calibrations”  needed  to  establish  the  various  composite  specific  properties 
and  damage  parameters  shown  in  Table  1.  The  effect  of  the  coating  in  coated  fibers  is 
included  in  the  transverse  modulus  of  the  fiber  listed  in  the  table.  The  table  also  lists  the 
references  in  which  the  test  data  analyzed  in  this  work  have  been  reported. 

Note  that  a  single  tensile  test  is  needed  to  determine  values  of  all  needed  constants 
relevant  to  the  laboratory  test  conditions.  Furthermore,  the  model  is  versatile  enough  to 
capture  the  significantly  different  deformation  characteristic  of  the  composites  with 
sufficient  accuracy.  While  the  tensile  stress-strain  curve  of  the  two  cross-ply 
Sylramic™-iBN/SiC  (silicon  carbide)  composites  shown  in  Figure-5  have  a  bilinear 
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profile,  the  cross-ply  Nicalon™/SiC  composite  shown  in  Figure-6  shows  nonlinear 
behavior  early  during  load-up.  The  tensile  response  of  the  latter  lacks  a  distinct  ‘knee’  in 
the  curve.  Figure-7  shows  a  second  stiffening  in  the  tensile  behavior  of  unidirectional 
Nicalon™/CAS  (calcium  alumino-silicate)  composite  close  to  failure.  The  model  does 
capture  this  behavior,  but  has  to  include  fiber  failure  to  accurately  model  the  observed 
data.  In  all  cases  it  was  found  that  the  median  strength  of  the  matrix,  <j0  ,  approximately 

coincided  with  the  stress  in  the  matrix  in  undamaged  composite  corresponding  to  the 
proportional  limit  in  the  stress-strain  curve  of  the  composite. 

Figure  8  shows  the  predicted  stress-strain  behavior  of  a  SiC/SiC  composite  when  it  is 
unloaded  from  certain  load  levels,  and  re-loaded.  Upon  unloading  the  composite  is 
predicted  to  have  linear  stress-strain  response  such  that  there  is  full  recovery  of  the  strain 
when  the  composite  is  fully  unloaded.  This  is  consistent  with  the  assumption  that  the 
nonlinearity  observed  during  tensile  loading  is  primarily  due  to  micromechanical  damage 
of  the  matrix.  The  prediction  using  the  above  model  required  no  additional  calibration 
beyond  that  described  earlier.  As  shown,  the  prediction  is  in  general  agreement  with 
measurements.  Unloading  portions  of  the  data  indicate  damage  accumulation  with 
increasing  load.  Similar  prediction  of  the  change  in  the  secant  modulus  of  unidirectional 
SiC/glass  composite  is  shown  in  Figure  9.  Again  there  is  good  agreement  between  the 
data  and  the  model  prediction. 

Figure  10(a)  shows  damage  state  at  room  temperature,  represented  by  the  parameters 
rj  and  p ,  plotted  against  increasing  applied  tension,  for  the  SiC/SiC  composite  009 
shown  in  Figure  5.  The  product  of  the  two  damage  parameters  is  also  plotted  in  the 
figure.  This  quantity,  which  is  the  product  of  the  transverse  and  longitudinal  normalized 
total  crack  lengths,  can  be  interpreted  as  a  measure  of  the  crack  density  in  the  matrix.  The 
plot  shows  that  below  a  critical  value  of  the  applied  stress  there  are  no  significant  cracks 
in  the  matrix.  Matrix  cracking  starts  as  soon  as  the  compressive  matrix  stress  in  the 
matrix  is  overcome  by  the  applied  stress.  The  value  of  the  applied  stress  when  this 
occurs  is  the  critical  ‘matrix  cracking  stress’  in  Figure  10(a).  Beyond  this  value  of  the 
load  the  micro  damage  in  the  matrix  progressively  increases  until  it  approaches 
saturation.  It  can  also  be  observed  in  Figure  10(a)  that  the  77  parameter  becomes  0.5  at 
applied  stress  of  166  MPa,  which  is  close  to  the  proportional  limit  of  the  009  composite 
in  Figure  5.  Thus  according  to  the  CMC  model  presented  here  the  density  of  the 
transverse  cracks  approaches  half  its  value  at  saturation  near  the  proportional  limit  of  the 
composite. 

The  product  of  the  damage  parameters  77  and  p  is  plotted  against  the  applied  stress  in 
Figure  10(b)  for  the  two  different  SiC/SiC  composites  shown  in  Figure  5.  This  quantity, 
which  is  the  product  of  the  transverse  and  longitudinal  normalized  total  crack  lengths,  can 
be  interpreted  as  a  measure  of  the  crack  density  in  the  matrix.  The  plot  also  shows  the 
linear  crack  density  estimated  from  the  acoustic  energy  (AE)  response  of  the  two 
composites  as  reported  in  reference  [19].  AE  is  a  measure  of  the  acoustic  signals 
associated  with  matrix  cracking  and  so  the  AE  curve  is  a  measure  of  the  number  of 
matrix  cracks  in  the  composite.  Thus  even  though  there  is  no  one-to-one  correlation 
between  the  quantity  77*77  and  the  AE  curve  they  can  be  compared  in  a  qualitative  sense 
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since  they  both  are  measures  of  the  growth  of  matrix  cracks  in  the  composite.  Figure 
10(b)  shows  that  the  AE  curve  and  the  product  of  77  and  p  both  have  a  sigmoidal  profile: 
they  both  start  slowly  and  increase  rapidly  before  approaching  a  saturation  value.  In  the 
CMC  model  cracking  starts  only  after  the  matrix  residual  stress  is  overcome.  Since  the 
damage  functions  p  and  77  of  the  above  CMC  model  are  nonnalized  their  product 
approaches  a  value  of  1.0  that  corresponds  to  fully  cracked  matrix.  Another  observation 
that  can  be  made  from  Figure  10(b)  is  that  the  quantity  p*t 7  preserves  the  material  ranking 
seen  in  the  AE  data.  Thus  the  matrix  in  the  009  material  is  predicted  by  the  model  to 
crack  more  rapidly  than  that  in  the  Oil  material  and  this  is  consistent  with  the  AE  data. 

The  growth  of  the  damage  parameters  77  and  p  in  the  Nicalon™/SiC  composite  along 
with  that  of  their  product  with  the  applied  stress  is  shown  in  Figure  11.  As  shown  in 
Figure  6  the  stress-strain  behavior  of  this  composite  exhibits  nonlinear  behavior  at  a  very 
small  load.  The  damage  functions  in  Figure  1 1  shows  that  the  damage  growth  starts  at  a 
much  lower  value  of  the  applied  stress  as  compared  to  the  SiC/SiC  composites  in  Figure 
10. 

Crack  density  in  SiC/glass  determined  by  counting  the  individual  cracks  in  the  specimen 
gage  area  are  reported  in  reference  [1].  Figure  12  compares  the  damage  state  prediction 
using  the  present  model,  represented  by  the  product  of  77  and  p ,  with  the  crack  density 
data  for  this  material.  In  this  case  too  the  model  captures  the  matrix  crack  growth 
behavior  observed  in  the  data. 

As  mentioned  earlier  with  the  exception  of  fiber  properties  and  architecture  specifics  of 
the  composite,  all  other  input  data  for  the  model  were  found  from  the  stress-strain  data. 
Due  to  the  highly  nonlinear  nature  of  the  problem  it  is  difficult  to  obtain  a  rigorous 
mathematical  proof  of  the  uniqueness  of  the  damage  constants  for  each  composite. 
However  we  can  use  the  closed-fonn  expression  (28)  determine  an  optimized  set  of 
damage  constants  that  best  agree  with  the  test  data.  All  the  damage  constants  reported  in 
this  work  were  detennined  in  this  manner. 

The  advantage  of  this  model  is  that  once  the  CMC  material  is  characterized  the  model 
can  be  used  to  study  the  effect  of  varying  geometry  and  loads.  Further  since  the  model  is 
described  in  tenns  of  the  basic  constituent  properties  it  can  also  be  used  to  estimate  the 
effect  of  changing  variables  such  as  porosity  and  fiber  volume  fraction.  The  predicted 
tensile  response  of  Sylramic™-iBN/SiC  CMCs  with  different  volume  fractions  is 
compared  with  test  data  from  reference  [19]  in  Figure  13.  The  model  predictions  were  all 
made  using  the  constituent  properties  determined  from  the  01 1-CMC  data  (Table- 1).  The 
CMCs  in  Figure  13  seem  to  differ  in  more  than  just  their  volume  fraction  as  indicated  by 
their  initial  modulus  not  changing  monotonically  with  the  fiber  content.  The  model  does 
however  capture  the  ranking  and  profile  of  the  overall  behavior  of  the  composites.  This 
application  of  the  above  CMC  model  illustrates  how  it  can  be  used  in  material  evaluation. 

Another  feature  of  the  model  is  the  coupling  of  the  constituent  material  inelasticity  with 
damage.  As  an  example  we  use  the  equations  to  model  the  time-dependent  behavior  of 
oxide/oxide  composites  in  which  the  matrix  cracks  easily  while  the  fiber  creeps  under 
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load.  Figure  14  shows  the  comparison  of  the  solution  from  the  model  with  test  data  [21] 
from  woven  cross-ply  composite  with  Nextel™720  fiber  and  alumina-silicate  (AS-0) 
matrix.  To  capture  the  observed  composite  behavior  we  assume  the  liber  strain  rate,  s‘f , 
to  follow  the  relation 

e‘f  =  1.6x1 0-12cr2f73t-OJ55  (31) 

Here  the  fiber  stress,  <jf,  is  expressed  in  MPa  and  time,  t,  is  in  seconds.  The  other 

material  properties,  including  damage  constants,  used  to  do  the  modeling  are  given  in 
Figure  14.  Thus  as  the  liber  creeps  and  the  matrix  cracks  they  shed  load  to  each  other  so 
that  both,  the  amount  of  creep  and  the  crack  density,  increase  continuously  with  time; 
though  in  this  example  the  increase  in  damage  strain  is  not  significant.  Figure  15  shows 
the  predicted  change  in  the  creep  and  damage  strains  up  to  3  hours  when  the  matrix 
becomes  fully  cracked  and  the  damage  strain  reaches  a  steady  state  value. 

Conclusion 

In  this  paper  we  present  a  mechanics  based  life-prediction  model  for  brittle  matrix 
composites  that  considers  the  deformation  due  to  dominant  damage  mechanisms  such  as 
matrix  microcracking  and  liber-matrix  debonding.  The  model  also  includes 
consideration  of  inelastic  deformation  of  one  of  the  constituents  such  as  creep  of  the 
liber.  All  the  material  properties  and  model  constants  can  be  obtained  by  charactering  a 
single  tension  and  single  creep  test  of  the  material  under  each  test  condition.  Once 
characterized  the  model  can  be  used  to  predict  the  behavior  under  other  geometry  and 
load  conditions.  The  model  also  detennines  the  damage  state  of  the  material  under  load. 
Examples  of  the  application  of  the  model  to  several  ceramic  matrix  composites  with 
similar  damage  mechanisms  have  been  presented.  This  modeling  approach  can  be  used 
in  as  a  basis  to  develop  a  comprehensive  life  prediction  model. 
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APPENDIX 


The  properties  of  the  effective  matrix  (m),  which  includes  the  actual  matrix  (M)  and  the 
transverse  fibers  in  a  balanced  cross-ply  composite,  are  given  by 


E.  = 


[E,-E,vf!2\ 


l-vf/'2 


a 


a...  = 


[Ef  vf/2  +  Em  (l-vf/2)\-Ef  a  f  vf/2 
Em(l~Vf/2) 


(Al) 

(A2) 


where  the  effective  properties  of  the  cross-ply  composite  along  a  fiber  direction  are  given 
by 


E  L  +  Et 
E,  =—± - 


CCj  = 


E  LaL  +  ETaT 

Ye, 


(A3) 

(A4) 


and  the  longitudinal  (L)  and  transverse  (T)  properties  of  a  single  ply  are  given  by 

El  =  Efvf+EM(l-vf)  (A5) 

(A6) 


EjjEm 

Et  =  t - - — -t - y 

EMVf+EA1-Vf). 


aL  = 


f  1  ^JTV  yf, 
EfafVf+EMaM(l-Vf) 


aT  = 


=  {l  +  vf  )afvf  +  {l  +  vM)aM(l-vf)-aL  [vfvf  +  vM  (l-vf  )J 


(A7) 

(A8) 
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Table- 1:  Properties  of  CMCs 


CMC 

1  HI 

120  [20] 

009  [19] 

011  [19] 

068  [19] 

Geometry 

Unidirectional 

0/90  PW1 

0/90  5H-SW2 

0/90  5H-SW2 

0/90  5H-SW2 

mmp 

0.35 

0.40 

0.00 

0.05 

HEM 

1000 

1450 

1450 

1450 

1450 

Fiber 

NicalonIM 

Nicalon"vl 

Sylramic"vl-BN 

SylramicIM-BN 

Syl  ramie  "VI-BN 

EfL  (GPa) 

190 

333 

333 

333 

EfT  (GPa) 

190 

180 

173 

173 

173 

Vf 

0.17 

0.17 

0.17 

0.17 

0.17 

Of  (p/° C) 

3.55 

3.90 

2.70 

2.70 

2.70 

Matrix 

CAS 

CVI-SiC 

Ml-SiC 

Ml-SiC 

Ml-SiC 

Em  (GPa) 

93 

295 

331 

181 

321 

vm 

0.125 

0.125 

0.125 

0.125 

0.125 

am  (M/°C) 

4.75 

3.50 

2.10 

2.10 

2.10 

c t0  (MPa) 

256 

47 

140 

158 

365 

In 

2.648 

0.410 

0.459 

0.518 

0.770 

Pc 

-0.119 

-0.050 

-0.117 

-0.190 

-1.674 

Pn 

3.945 

1.200 

1.602 

1.746 

2.867 

1  PW  -  Plain  Weave 

2  5H-SW  -  5  Harness  Satin  Weave 
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tr 


Load  Direction 


1  -  Minimal  or  no  additional  damage 

2  -  Intra-tow  cracking  in  transverse  tows 

3  -  Inter-tow  cracking,  some  splitting 

4  -  Splitting,  microcrack  growth 


Figure  1(a):  Response  of  SiC/SiC  CMC  under  uniaxial  tension. 


Time 


Fiber  creep  &lnitial  matrix 
cracking  in  transverse  tows 

Environmentally  assisted 
progressive  matrix  cracking. 

Splitting/Interfacial  sliding, 
followed  by  fiber  fracture 


Figure  1(b):  Response  of  SiC/SiC  CMC  under  sustained  loading. 
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Fiber  Material: 

Elastic/Inelastic 

Matrix  Material: 

Elastic 


a=Debond  Length 
s=Shear  Lag  Length 


a  +  s 


_  Number  of  Type  B  Cells 
^  Number  of  Type  U  Cells 


® 

J2 

Matrix 

a) 

-Q 

LL 

LL 

Region  B, 

1 

s 

]  Region  B2 

1_ 

a) 

■Q 

U  L. 

Matrix  £ 

Ll 

LL 

1 

a 

T 


Figure  2:  Unit  Cells  in  CMC  Damage  Model. 
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Composite) 


Figure  3:  CMC  Damage  Model. 
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Figure  4(a):  Comparison  of  Numerical  Solution  with  Exact  Solution  of  CMC  Model. 


0  100  200  300  400 

Stress  (MPa) 


Figure  4(b):  Evolution  of  Damage  Functions  in  CMC  Model. 
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Strain  (m/m) 


Figure  5:  Tensile  Response  of  Sylramic-iBN/MI-SiC  Cross-Ply  CMCs. 


Strain  (m/m) 


Figure  6:  Tensile  Response  of  Nicalon/CVI-SiC  Cross-Ply  CMC. 
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Strain  (m/m) 


Figure  7:  Tensile  Response  of  Nicalon/CAS  Unidirectional  CMC. 
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Figure  8:  Tensile  Load-Unload  Response  of  Sylramic-iBN/MI-SiC  Cross-Ply  CMC. 
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Strain  (m/m) 


Figure  9:  Change  in  Secant  Modulus  of  Nicalon/CAS  Unidirectional  CMC. 


Figure  10(a):  Damage  Evolution  in  Sylramic-iBN/MI-SiC  Cross-Ply  CMC. 
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Figure  10(b):  Comparison  of  Damage  Evolution  in  Sylramic/MI-SiC  Cross-Ply  CMCs. 
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Figure  11:  Damage  Evolution  in  Nicalon/CVI-SiC  Cross-Ply  CMC. 
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Composite  Stress  (MPa) 


Figure  12:  Damage  Evolution  in  Nicalon/CAS  Unidirectional  CMC. 


Strain  (m/m) 


Figure  13:  Effect  of  Volume  Fraction  on  Sylramic-iBN/MI-SiC  Cross-Ply  CMC. 
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Figure  14:  Creep  Deformation  in  Nextel/AS-0  Cross-Ply  CMC. 


Figure  15:  Change  in  Strain  Components  During  Creep  of  Nextel/AS-0  Cross-Ply  CMC. 
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